Random unitary maps for quantum state reconstruction 
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We study the possibility of performing quantum state reconstruction from a measurement record 
that is obtained as a sequence of expectation values of a Hermitian operator evolving under repeated 
application of a single random unitary map, Uq. We show that while this single-parameter orbit 
in operator space is not informationally complete, it can be used to yield surprisingly high-fidelity 
reconstruction. For a d-dimensional Hilbert space with the initial observable in su(d), the mea- 
surement record lacks information about a matrix subspace of dimension > d — 2 out of the total 
dimension d 2 — 1. We determine the conditions on Uo such that the bound is saturated, and show 
they are achieved by almost all pseudorandom unitary matrices. When we further impose the con- 
straint that the physical density matrix must be positive, we obtain even higher fidelity than that 
predicted from the missing subspace. With prior knowledge that the state is pure, the reconstruction 
will be perfect (in the limit of vanishing noise) and for arbitrary mixed states, the fidelity is over 
0.96, even for small d, and reaching T > 0.99 for d > 9. We also study the implementation of this 
protocol based on the relationship between random matrices and quantum chaos. We show that the 
Floquet operator of the quantum kicked top provides a means of generating the required type of 
measurement record, with implications on the relationship between quantum chaos and information 
gain. 

PACS numbers: 32.80.Qk,42.50.-p,02.30.Yy 



I. INTRODUCTION 



Quantum state reconstruction (QSR, or quantum to- 
mography) is a fundamental tool in quantum information 
science that has been carried out in a variety of systems 
and in a variety of protocols [1, 2]. The essential proce- 
dure of QSR is to use the statistics of the measurement 
results on an ensemble of identical systems to make a best 
estimate of the prepared state po- This can be achieved, 
e.g., through a series of strong projective measurements 
of a set of Hermitian observables [3] or through contin- 
uous weak measurement of a time-series of observables 
[4, >]. High fidelity QSR typically requires an "informa- 
tionally complete" measurement record. One can obtain 
informational completeness by measuring the expecta- 
tion values of a set of Hermitian operators that span 
an operator basis for po, or some more general opera- 
tor "frame" [6, 7]. Restricting our attention to a Hilbert 
space of finite dimension d, and fixing the normalization 
of po, the set of Hermitian operators must form a ba- 
sis for the Lie algebra su(d). Laboratory realization of 
such a record is intimately tied to controllability, i.e., the 
ability to reconfigure the apparatus in such a way as to 
generate arbitrary unitary maps. In the continuous mea- 
surement context, when the system is controllable it is 
possible to choose control fields for the system such that, 
when viewed in the Heisenberg picture, the observables 
evolve over the span of the algebra. While the necessity of 
information completeness is rigorous if one requires high 
fidelity for the reconstruction of all arbitrary states, in a 
number of situations this condition can be substantially 



relaxed. Examples include schemes that are designed to 
achieve high performance of the reconstruction on aver- 
age [8] , or over only some restricted set of the state space 
[')]. For these protocols, the performance of a restricted 
set of measurements is often nearly as good as for an 
informationally complete set of measurements, and yet 
require dramatically fewer measurement resources. 

In this paper, we study another example of informa- 
tionally incomplete measurements that nonetheless can 
be used in a high-fidelity QSR — measurement of a time- 
series of operators generated by a single-parameter ran- 
dom evolution. Consider continuous weak measurement 
of an observable, O, through a meter that couples to an 
ensemble of N identical systems. The members of the en- 
semble undergo identical, separable time evolution in a 
well chosen manner. Assuming the subsytems remain in a 
product state, or equivalcntly that a mean-field approx- 
imation is valid, we can write our measurement record 
quite generally as 



M(t)=N(0) (t) + SM(t), 



(1) 



where, 8M(t) describes the deviation from the mean 
value arising from state-dependent quantum uncertainty 
and noise in the detection system. For unitary evolution, 
(O) (t) = Tr (W(t)OU(t)p ). The QSR problem is to 
retrodict the initial state of the system po from the sig- 
nal M. (t) . We can simplify the analysis of this problem 
by considering a discrete set of measurements at intervals 
At, {O n = W(nAt)OU(nAt)}. The ultimate fidelity of 
the QSR will be limited by the finite signal-to-noise ratio. 
While the choice of unitary evolution necessary to deter- 
mine an arbitrary po from M. n is not unique, a necessary 
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and sufficient condition is that the set {£>„} be infor- 
mationally complete. A good strategy is to choose the 
dynamics such that for each n, U (nAt) is a random ma- 
trix, chosen from an appropriate Haar measure. In that 
case our measurement record is not only provably infor- 
mationally complete but is also unbiased over time. Sup- 
pose, however, we choose U{nAt) = (Uo) n , where Uq is 
a fixed unitary matrix. In this case, the observable series 
O n traces out a single orbit in operator space; we call this 
a one-parameter measurement record. As we will show, 
the record is not informationally complete, but neverthe- 
less can lead to high fidelity QSR for all states but a set 
of small measure if Uq is a random unitary, especially 
for large dimensional spaces. These results elucidate the 
connection between random evolution and information 
gain at the quantum level. 

The remainder of the article is organized as follows. 
In Sec. II we show that the measurement operators gen- 
erated from a single parameter trajectory cannot span 
the entirety of the operator algebra, su(d) , but that the 
operators that lie outside of the subspace of the mea- 
surement record are a vanishingly small fraction in the 
limit d — ► oo. Next, we study the performance of QSR 
using the weak measurement protocol [4] for these incom- 
plete measurement records. We show that even at small 
d, when one includes the physical constraint of positivity 
of the density matrix, QSR performs surprising well for 
almost all quantum states, well beyond that expected if 
one had solely considered the vector space geometry of 
the Lie algebra. Finally, in Sec. Ill, we connect these 
abstract results to physical realizations using the unitary 
Floquct maps of the quantum kicked top [10] whose asso- 
ciated classical dynamics is chaotic. As quantum chaos 
is associated with pseudorandom matrix statistics, this 
protocol provides intriguing new signatures of quantum 
chaos in QSR. 



II. ONE-PARAMETER MEASUREMENT 
RECORDS 



In this section, we study whether or not informa- 
tion completeness for QSR is achievable from a one- 
parameter measurement record. The one-parameter or- 
bit in operator space is defined by the time-series O n = 
{Ul) n O{U a ) n : where O is a Hermitian operator and Uq 
is a fixed unitary matrix. We will restrict the observable 
O to have zero trace since the component proportional 
to the identity gives no useful information in QSR. We 
thus ask, is it possible to reconstruct a generic quan- 
tum state po if one can measure the expectation values 
of all of the observablcs in the time series? To answer 
this, we consider A = span{0„}, and determine the size 
the orthocomplement subspace with respect to the trace 
inner product, A±; operators in this set are not mea- 
sured in the time-series. Such missing information ren- 
ders the measurement incomplete, and thus incompatible 
with perfect QSR, no matter what signal-to-noise ratio 



is available in the laboratory. 

To find the dimension of A, consider the sub- 
space of operators that are preserved under conjuga- 
tion by Uq, Q = | .9 G 5u(d) U cjUq — . Let B = 

{g EG\ Tr(gO) = 0} . It thus follows that B C A± since 
VgeB 

Tr(O n g) = Tr {(U^) n O(U ) n g) = Tt{Og) = 0. (2) 

As the two spaces are orthogonal, dim A + dim ,8 < 
dim(su(e?)) = d 2 — 1. Now, if Uq has nondegener- 
ate eigenvalues, Q will be isomorphic to the the largest 
commuting subalgebra of su(d) (the Cartan subalgebra), 
but for degenerate Uq, Q will contain additional ele- 
ments. Since the Cartan subalgebra has dimension d — 1, 
dim Q > d — 1. By definition, B is obtained from Q by 
projecting out one direction in operator space, and thus 
dim ,8 = dimC/ — 1 > d — 2. It follows that 

dim .4 < dim(su(d)) - dim 6 < d 2 - d + 1. (3) 

This is the first principal result - a one-parameter mea- 
surement record is not informationally complete since 
dim.A_L > (when d > 2). But, it remains to be seen 
how much the missing information impacts the fidelity of 
QSR. An immediate question is to determine the condi- 
tions on Uq and O required to saturate bound in Eq. (3). 
Since Uq is a unitary matrix it is always diagonalizablc 
as 



c/b = |i>0l, 

i=i 

and in this basis O n has the representation 



(4) 



°n= £ e- m ^-^(k\0\j)\k)(j\. (5) 
j,k=i 

The diagonal component has no 7i-dependence, so it is 
useful to rewrite O n as 

d d 

On = £010 li) li) 01 + £ e^-M ( k \ O \j) \k) oi . 

3=1 j^k 

(6) 

To show that that A is spanned by d 2 — d + 1 linearly 
independent matrices we must have that 

d 2 ~d 

£ a n O n = iff a n = Vn. (7) 

n=0 

We can write this condition out explicitly using Eq. (6) 



as 



'd 2 -d \ d 

£a n b*>0'I + 

n=Q J j = l 

E[E«-e^-«) (k\O\j)\k){j\ = 0. (8) 

jjtk \n=0 I 
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The system of equations is underconstrained if cither 
U\0\j) = for all j or (j\0\k) = for any j ^ k. 
Assuming this is not the case, the condition for linear 
dependence is given by a set of linear equations on a„ of 
the form 



(l 



1 x d 2 -d x 



2 

d 2 -d 



-d 



-d 

-d / 



a 
a i 

a d 2 -d 



0. (9) 



Here we have written xq = 1 and x m = e~ t ^ i ~^ k ' , for 
some indexing of the pairs (j, k) to 1 < m < d 2 — d. The 
condition for linear independence is simply detV =/= 0. 
Expressed as above, one can see that V is an instance 
of a Vandermonde matrix, whose determinant is easy to 
evaluate through the formula [11] 



0<j<k<d 2 -d 



(10) 



For our system of equations to become linearly dependent 
before saturating the previous bound, we would need that 
e -i(4>j-(p k ) _ e -i(<t>j> -<t>k>) f or some distinct pair of the 



or e 







1 for some 



couples (j, k) and 
U,k). 

In summary, in order for the dimension of the span of a 
one-parmeter measurement record to saturate the bound 
of dim .4 = d 2 — d+1, the eigenphases, <pj, and the eigen- 
vectors, \j), of U must satisfy the following constraints: 



1. 3 j s.t., 

2. Vj^fc, 
3- Vj^j', 



(k\ O \j) ± 

<t>i - <t>k ^ (j>j< 



(11) 



Note that the third condition enforces that both the 
eigenphases, as well as their pair wise differences, must 
be distinct. There is an interesting interpretation of 
these conditions from the perspective of universal con- 
trol. While f/o defines a one-parameter trajectory, the 
set {Uo.e 10 } defines a universal set of unitary matrices 
that can generate arbitrary maps, if the conditions in 
Eq. (11) are satisfied [12]. 

If Uq is a unitary matrix chosen randomly from the 
Haar measure on SU(d), the saturation conditions will 
almost surely be satisfied, independent of O. Therefore, 
a generic unitary evolution will almost always generate a 
measurement record that spans the full d 2 — d + 1 op- 
erators. In fact, the typical members of many types 
of pseudorandom ensembles of unitary matrices satisfy 
these constraints, e.g. unitaries drawn from f-designs 
or approximate i-designs [13, 14, 15, 16], random quan- 
tum circuits [17, 18], as well as unitary evolutions that 
possess globally chaotic dynamics in the classical limit 
[10, 19, 20]. These types of pseudorandom evolutions are 
more readily available in practical situations, providing 



possible avenues to test these results in laboratory im- 
plementations. 

The results of this section show that a one-parameter 
evolution generates a measurement record that misses a 
subspace of dimension d — 2 out of the full su(d) algebra 
whose dimension is d 2 — 1. For very large Hilbert space 
dimensions, the implication is that all but a vanishing 
fraction of the information regarding measurements on 
the quantum system is contained in this type of record. 
The reconstruction fidelity is a fairly complicated nonlin- 
ear function of the measurement record, 




(12) 



It is not clear that this fidelity will be directly related 
to the fraction of operator space spanned by our set of 
observables. Surprisingly, the situation is in fact more 
favorable than this naive assumption. In the next sec- 
tion we will see that merely requiring the reconstructed 
density matrix be positive provides a powerful constraint, 
allowing us to use a one-parameter measurement record 
induced by the orbit of a single pseudorandom unitary 
matrix to perform very high fidelity reconstructions even 
for small dimensional Hilbert spaces, for all but a very 
small subset of states. 



III. DENSITY MATRIX RECONSTRUCTION 
FROM AN INCOMPLETE MEASUREMENT 

The QSR protocol we consider was first proposed by 
Silberfarb et al. [4] and implemented by Smith et al. 
[5]. In this scenario one has access to an ensemble of N 
identically prepared systems all initialized to the same 
state, po- The system is weakly measured yielding the 
record given in Eq. (1). For sufficiently weak coupling, 
the deviation of the measurement result from the quan- 
tum expectation value is dominated by the noise on the 
detector (e.g., shot noise of a laser probe) rather than 
the quantum fluctuations of measurement outcomes in- 
trinsic to the state (known as projection noise). In this 
case, there is negligible backaction on the quantum state 
during the course of the measurement, and the ensemble 
remains factorized. We treat the detector noise as Gaus- 
sian white noise defined as 5M(t) = crW(t), where W(t) 
is a Weiner process and a determines the noise variance. 

We examine a stroboscopic time-series of the measure- 
ment record, where the observables evolve according to 
the one-parameter trajectory discussed in the previous 
section. At discrete times t = nAt, 



M n = NTr(O nPo )+aW n , 



(13) 



where O n and po are Heisenberg operators. The prob- 
lem is thus reduced to one of stochastic estimation of 
po given {M n }. To accomplish this, the density op- 
erator is expanded in a basis of Hcrmitian matrices, 
{E a }, so po = J2 a r aE a + I/d. We take po to be unit 
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trace, explicitly removing the identity from the opera- 
tor basis. Defining a rectangular matrix with elements 
O na = Ti(O n E a ), the measurement time-series given in 
Eq. (13) can be expressed as the vector M = NOr+aW . 
Because the fluctuations around the mean are Gaussian 
distributed, the maximum-likelihood-estimate of the un- 
known parameters {r a } is the least-squared fit, given by 
the vector 

r ML = }_(& r d)- 1 6 T M. (14) 

If the measurement record is informationally complete, 
which occurs when the covariance matrix C = 6 T 6 has 
full rank d — 1, and in the absence of measurement noise, 
the maximum-likely-estimate, p ML = ^2 a r^ L E a + I/d, 
is exactly po- If the set of observables {O n } is not infor- 
mationally complete, then C is not full rank and we must 
replace the inverse in Eq. (14) with the Moore-Penrose 
pseudo-inverse, i.e., inverting only over the space in which 
the covariance matrix has support. In this case, a state 
with support on the null space of C will yield a maximum 
likelihood estimate with sub-unit reconstruction fidelity. 

In the presence of measurement noise, or when the 
measurement record is incomplete, the estimate p ML can 
have negative eigenvalues which are unphysical. To ob- 
tain a better reconstruction of the state, we seek the 
physical density matrix "closest" to p ML . We use the 
covariance matrix C as a cost function or metric to mea- 
sure the distance between p ML and new estimate p by 
defining 

|| r ML _ F ||2 = ( r ML _ f )T c(r ML _ f) ^ 

Technically speaking, this quantity is not a norm but 
rather a scminorm, meaning that there exist some vectors 
v such that j|v|| = but v 7^ 0. This cost penalizes 
us for taking displacements in directions that increase 
the variance in the measurement uncertainty. We thus 
need an f that minimizes the distance ||r ML — f || while 
respecting the constraint ^2i a f a E a + I/d > 0. While 
there is generally no analytic solution to this problem, the 
optimization is a semidefinite program which is efficiently 
solvable numerically [21, 22]. 

The condition of positivity is a powerful constraint that 
describes correlations between observables that can lie 
along orthogonal directions in operator space. For exam- 
ple, in the case of a 2-level quantum system, if (a z ) = 1, 
positivity implies (a x ) = (a y ) = 0, fully specifying the 
state from a single expectation value. In the context of 
noisy measurcmets, the positivity constraint allows us 
to perform high fidelity QSR in the face of uncertainty 
by enforcing consistency conditions on our measurement 
outcomes. When we consider incomplete measurement 
records, positivity can place bounds on the means of ob- 
servables which otherwise would be completely undeter- 
mined. This can greatly increase the fidelity of QSR. In- 
tuitively, while many vectors r might minimize Eq. (15), 
only very few of these arc also compatible with positivity. 
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FIG. 1: (Color online) Numerical simulations the QSR pro- 
tocol for pure states as function of n th expectation value 
measured in the time-series, and for different dimensions of 
the Hilbert space, d. Each data point represents the aver- 
age reconstruction fidelity of 100 pure states drawn from the 
Fubini-Study measure, additionally averaged over measure- 
ment records derived from ten different Harr-random unitary 
propagators. 



As we show below, in the context of a one-parameter mea- 
surement record generated by a single random matrix, 
the requirement of positivity provides substantial addi- 
tional information leading to very high fidelity of QSR, 
well beyond what one would naively predict. 

We are now prepared to quantitatively analyze the per- 
formance of our QSR protocol in the case of the one- 
parameter measurement record arising from an incom- 
plete set of observables that satisfy Eq. (11). As our 
system, we consider a system with spin J described by 
a Hilbert space of dimension d = 2 J + 1. We will fix 
O = J z and select a random unitary matrix Uq from the 
Haar measure on SU(d). Such a random matrix will al- 
most always satisfy the constraints of Eq. (11), except for 
a set of measure zero. As our goal is to determine how the 
information missing in a subspace of observables impacts 
the QSR fidelity, we will simplify the analysis by assum- 
ing that noise on the measurement is vanishingly small. 
We study the performance of different classes of states, 
randomly chosen by an appropriate measure. For each 
set of states, we will look at the average fidelity between 
the initial and reconstructed states, (J 7 ) = J dpoJ-(p, po), 
where dpo is a measure on the space of density operators. 

The simplest case to analyze is when we have prior in- 
formation that po is a pure state, li/'o) ■ Figure 1 shows the 
average fidelity as one sequentially measures the expecta- 
tion value of the n th observable in the series, for different 
dimensions of the Hilbert space d. Averages are taken 
for 10 choices of random unitary matrices, each of which 
is averaged over 100 random pure states distributed on 
the Fubini-Study measure [23]. Two striking features are 
seen in these plots: (i) unit fidelity is achieved for any 
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d even though the record was said to be informationally 
incomplete; (ii) the protocol reconstructs the state well 
before we measure all d 2 — d + 1 independent observ- 
ables. The inclusion of positivity dramatically improves 
the reconstruction fidelity for pure states. In fact, a one- 
parameter measurement record generated by a random 
Uq can be used to reconstruct almost all pure states per- 
fectly in the absence of noise. 

The performance of the QSR protocol can be under- 
stood given the prior information we have assumed. A 
pure state is specified by 2d — 2 real parameters, whereas 
we measure d 2 — d+ 1 expectation values. Thus, it should 
come as no surprise that the measurement record con- 
tains enough information to reconstruct the state. In 
fact, in this case one can use positivity to explicitly re- 
cover the missing information exactly from the measure- 
ment record, without resorting to the numerical semidcf- 
inite program discussed above. In general, the missing 
information is associated with matrices that commute 
with Uq. Thus, when expressed in the eigenbasis of Uo, 
only the diagonal matrix elements of the density opera- 
tor might not be estimated. A necessary (but not gen- 
erally sufficient) condition for a matrix po to be positive 
semidefinite is that its matrix elements must satisfy the 
following set of inequalities: pupjj — \pij\ 2 > 0, i.e. all 
of the 2x2 matrix minors must be positive semidef- 
inite [11]. If additionally the state is pure, these in- 
equalities become equalities. Therefore, if any of the off- 
diagonal matrix elements are nonzero, we can completely 
determine all of the diagonal elements via the equations 
Pa = {\Pij\\Pik\) /\Pjk\- A special case is if all of the off- 
diagonal matrix elements of po are zero. Then the state 
must be one of the eigenvectors of Uo, but such states 
lie in a set of measure zero. Measurements that are in- 
formationally complete solely for pure states are called 
PSI-compete [24, 25]. 

While one can easily explain the high-fidelity perfor- 
mance of the QSR protocol in the case of a pure state, for 
mixed states, this is far from clear, and the power of the 
positivity constraint comes fully to the fore. For mixed 
states, and d > 2, the average fidelity of our reconstruc- 
tion will never reach unity because some density matri- 
ces cannot be reconstructed from the information in our 
incomplete measurement record. For example, some con- 
vex combinations of eigenstates of Uq are indistinguish- 
able from the maximally mixed state, even though the 
fidelity between the two can be very small. Nonetheless, 
as we see below, the one-parameter measurement record 
generated by a single random unitary still performs very 
well on average, even for generic mixed states. 

Figure 2 shows the average fidelity for two choices of 
measures on density matrices, the Bures measure and the 
Hilbert-Schmidt measure [23, 26], with states sampled 
according to the construction provided by Osipov, Som- 
mers and Zyczkowski [27]. For both distributions we look 
at a long-time limit of the time-series, here 10(g? 2 — d+1) 
measurement steps, and plot the average fidelity as a 
function of the dimension of the Hubert space rather 
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FIG. 2: (Color online) Fidelity of the QSR protocol for mixed 
states as a function of the dimension of Hilbert space d. Each 
point represents the average reconstruction fidelity from a 
measurement record of length 10(d 2 —d+1) (a long-time 
limit) generated from a different random unitary propagator 
from the Harr measure. We average over 200 density matri- 
ces drawn from the Bures (blue crosses) or Hilbert-Schmidt 
(green x's) measures. For each dimension we show the average 
fidelities from twenty of these measurement procedures. 



than n. In the limiting case of negligible noise on the 
measurement, we have already extracted all the possible 
information about the state after d 2 — d+1 measurements. 
In practice, increasing the measurement record serves to 
smear out the information over the measured observablcs, 
leading to a more uniform distribution for the non-zero 
eigenvalues of C, which is numerically favorable. As seen 
in these two plots, on average, the one-parameter mea- 
surement records perform surprisingly well. In all cases 
the mean fidelity is greater than 0.96 with a minimum 
around d = 3 or d = 4. After this dip, the minimum of 
the fidelity looks to be monotonically increasing with the 
size of the Hilbert space. Additionally, the particular in- 
stantiation of the random unitary map appears to make 
very little difference (less than 0.01 fidelity), with the 
residual difference decreasing as the dimension increases. 



IV. EXAMPLE: QUANTUM KICKED-TOP 

In Sec. II we discussed that the conditions given in 
Eq. (11) can be satisfied by a pseudorandom unitary ma- 
trix, instead of a true random matrix sampled from the 
Haar measure. One such class of matrices are the Flo- 
quet maps associated with "quantum chaos", i.e., peri- 
odic maps whose classical dynamical description shows a 
globally chaotic phase space. An example is the quantum 
kicked-top [10], a system that recently has been realized 
in a cold atomic ensemble [28]. In this section we explore 
how our QSR protocol performs in this context, providing 
a possible route to laboratory studies, and novel signa- 
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tures of chaos in quantum information. 

The standard quantum kicked top (QKT) dynamics 
consists of a constant quadratic twisting of a spin ( "top" ) , 
punctuated by a periodic train of delta-kicks of the spin 
around an orthogonal axis. The Floquet operator for 
this perodic map is typically written as the product of 
noncommuting unitary matrices, 

t/ Q KT = e-^ 2 / J e-^. (16) 

The parameters 9 and <fi represent the angles of linear and 
nonlinear rotation respectively The dynamics exhibit a 
classically chaotic phase space for an appropriate choice 
of these parameters [10]. The connection between chaos 
in this system and random matrices has been well stud- 
ied, particularly, the relationship between the level statis- 
tics of the Floquet eigenvalues, chaos, and symmetry. 
Floquet maps associated with global chaos are random 
matrices that divide into different classes. If the Floquet 
operator is time-reversal invariant, the level statistics are 
that of the circular orthogonal ensemble (COE); with- 
out additional symmetry they are members of the cir- 
cular unitary ensemble (CUE). The latter group is U(d) 
or SU(cZ) depending on the context. The measurement 
records generated from matrices chosen from either the 
COE or CUE will satisfy the eigenvalue conditions de- 
scribed in Eq. (11) almost surely. The QKT is known 
to have a time-reversal symmetry and classically chaotic 
dynamics. It still does not have COE statistics in the full 
d = 2 J + 1 Hilbert space, however, due to an additional 
symmetry The QKT map is invariant under a 7r-rotation 
about the x-axis, leading to a parity symmetry. This sys- 
tem therefore has a doubly degenerate eigenspectrum, 
breaking the conditions in Eq. (11). While such Flo- 
quet operators generate a measurement record that have 
much less information relative to an arbitrary state, we 
can perform high-fidelity QSR in for states restricted to 
a subspacc defined by the additional symmetry, here, the 
states that have even parity under reflection around the 
x-axis. To do this we require that our initial operator O 
is also symmetric under reflection, e.g., O = J x . 

We present examples of this type of reconstruction in 
Fig. 3. Here we look at the QKT dynamics for a spin 
J = 3 particle (a d = 7 dimensional Hilbert space). We 
choose the parameters <j> = 7 and 9 = 0.228, values for 
which the classical phase space is known to be globally 
chaotic. Additionally, we let the noise on the measure- 
ment approach zero. For general mixed and pure states 
sampled from the full Hilbert space, this reconstruction 
performs poorly. However, if we restrict our attention 
to states that are eigenstates of parity, we reconstruct 
with near unit fidelity. As a check, we find that covari- 
ance matrix C has rank 19. This agrees with our predic- 
tions, since this space has a 4-fold degenerate — 1 parity- 
eigenspace and a 3-fold degenerate +1 parity-eigenspace. 
The measurement operators that preserve this symmetry 
must be block diagonal, lacking the (3 x 4) components 
from each of the two off-diagonal blocks. For this dimen- 
sion d 2 - d + 1 - 24 = 19. 
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FIG. 3: (Color online) Reconstruction fidelity for a variety 
of states versus n. The initial measurement observable is J x , 
and each subsequent observable whose expectation value we 
measure is obtained by evolving under the Floquet map of a 
quantum kicked top, Eq. (16). For generic random states in 
the whole Hilbert space (pure or mixed), the reconstruction 
preforms poorly. Density matrices that are invariant with re- 
spect to 7r-rotation around the x-axis, such as the cat state, 
an eigenstate of the parity operator, or an incoherent mix- 
ture of odd-parity J x eigenstates, are reconstructed with high 
fidelity. 
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FIG. 4: (Color online) Same as Fig. 3 when the evolution is 
given by the double kicked top. All states, pure or mixed, 
asymptote to fidelities near unity, with the pure states reach- 
ing their maxima quicker than the mixed states. 



We can examine the effects of a pseudorandom unitary 
that satisfies Eq. (11) on the whole space if we look at the 
"double kicked top" where we alternate kicking about x 
and y. Here the Floquet operator has the form 

U 2 kt = e-^l^e-^-e-^'^^e-^^. (17) 

If we choose = 0' = 6, 6 X = tt/2 and 9 y = 0.228, 
these operators have no time-reversal symmetry, or in- 
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deed any other symmetry, and so approximate the spec- 
trum of the CUE. Since the double kicked top Floquet 
operator shares no symmetries with J z wc can choose 
O = J z in order to satisfy the first condition of Eq. (11). 

In Fig. 4 we show the QSR performance for a time- 
series generated by the double kicked top in a Hilbcrt 
space of a spin-3 particle. Here the fidelity asymptotes 
to unity for all of our choices of states, as we would ex- 
pect from the simulations in Sec. III. The reconstruc- 
tion reaches its asymptotic fidelity when we have made 
d 2 — d + 1 = 43 measurements. As we saw previously, 
the pure states require less measurements to reach their 
asymptotic value and the QSR can be perfect in the ab- 
sence of noise. 



V. SUMMARY AND OUTLOOK 

We have studied measurement records that are derived 
by stroboscopically measuring the expectation values of 
a single observable of a system that is evolving under the 
repeated application of a single unitary map. We have 
shown that this record never contains complete informa- 
tion about the quantum state, however for unitary maps 
chosen randomly or pseudorandomly, only a vanishing 
fraction of the information is missing. When combined 
with the constraint of positivity, this incomplete mea- 
surement record led to a protocol for quantum state re- 
construction that had high-fidelity performance for typ- 
ical mixed and pure quantum states. For pure states we 
can achieve unit fidelity reconstruction (in the absence 
of noise) and for mixed states the fidelity is greater than 
0.99 for d > 9. 

A particular set of pseudorandom matrices we stud- 
ied in some detail are the Floquet maps generated by 
the quantum delta-kicked top. In the general case, these 
maps appear to be equally as effective for reconstruction 
as Harr-random unitary maps. In cases where the kicked 
top map exhibited additional symmetry, we were able 
to see that our reconstruction protocol required the ex- 
tra constraint that the measurement operator and states 
shared the symmetry as well. A map that is chaotic on 
the whole phase space saturates the bound of pseudo- 
random unitary operators on the whole Hilbert space. 



It is surprising that such a simple measurement pro- 
tocol should lead to such good average reconstruction 
fidelities. The reason for this appears to be a combina- 
tion of the mixing power of random evolutions and the 
constraints on state space associated with positivity. We 
do not yet have a rigorous explanation of these results, 
however, because the set of positive operators is a con- 
vex cone rather than a vector subspace [23], and thus 
it is difficult to quantify the volume of states that have 
both large support in the missing subspace of su(d) and 
are positive. Our conjecture is that both the increas- 
ing average fidelity and the decreasing dependence on 
the sampled unitary can be a explained based on con- 
centration of measure in analogy with the work in [9]. 
Essentially, most randomly sampled states have very lit- 
tle support on a the subspace that we do not measure. 
Related work in the field of matrix completion has shown 
that such "incoherence" between states and (incomplete) 
measurements can provably lead to high-fidelity state re- 
constructions especially when the states in question are 
low rank [29] or near to low-rank states [9]. Irrespective 
of a rigorous proof, it is an empirical fact that our pro- 
tocol works well with typical states and typical unitary 
maps. We expect that as the dimension of the Hilbcrt 
space increases, almost all of the states and unitary evo- 
lutions that are sampled will be very close to a typical 
value, resulting in high QSR fidelity. 



Acknowledgments 

We thank Andrew Scott for helpful discussions. STM 
was supported by Industry Canada and Quantum Works. 
STF was supported by the Perimeter Institute for The- 
oretical Physics. Research at Perimeter is supported by 
the Government of Canada through Industry Canada and 
by the Province of Ontario through the Ministry of Re- 
search & Innovation. STF thanks the Kavli Institute for 
Theoretical Physics where this work was completed; this 
research is also supported in part by the National Sci- 
ence Foundation under Grant No. PHY05-51164. IHD 
and CAR where support by the National Science Foun- 
dation, Grant 0903692. 



[1] M. G. A. Paris and J. Rehaek, Quantum State Estimation 
(Springer, 2004). 

[2] G. M. DAriano, M. G. A. Paris, and M. F. Sacchi, Ad- 
vances in Imaging and Electron Physics 128, 205 (2003). 

[3] G. Klose, G. Smith, and P. S. Jessen, Phys. Rev. Lett. 
86, 4721 (2001). 

[4] A. Silberfarb, P. Jessen, and I. Deutsch, Phys. Rev. Lett. 
95, 030402 (2005). 

[5] G. Smith, A. Silberfarb, I. Deutsch, and P. Jessen, Phys. 
Rev. Lett. 97, 180403 (2006). 

[6] G. M. DAriano, L. Maccone, and M. G. A. Paris, 



J. Phys. A: Math. Gen. 34, 93 (2001). 

[7] A. Roy and A. J. Scott, J. Math. Phys. 48, 072110 (2007). 

[8] S. Aaronson, Proc. Roy. Soc. A 463, 2088 (2007). 

[9] D. Gross, Y.-K. Liu, S. T. Flammia, S. Becker, and J. 
Eisert, e-print arXiv:0909.3304 (2009). 
[10] F. Haake, Quantum Signatures of Chaos (Springer- 

Verlag, Berlin, 2006). 
[11] R. A. Horn and C. R. Johnson, Matrix Analysis (Cam- 
bridge University Press, 1985), see Sec. 0.9.11 for Vander- 
monde determinants and Thm. 7.2.5 for positive semidef- 
inite criteria. 



8 



[12] C. Altafini, J. Math. Phys. 43, 2051 (2002). 

[13] A. Ambainis and J. Emerson, in Proc. of the 22nd Annual 

IEEE Conf. on Computational Complexity, pp. 129-140 

(2007). 

[14] D. Gross, K. Audenaert, and J. Eisert, J. Math. Phys. 

48, 052104 (2007). 
[15] C. Dankert, R. Cleve, J. Emerson, and E. Livine, Phys. 

Rev. A 80, 012304 (2009). 
[16] A. Roy and A. Scott, Designs, Codes and Cryptography 

53, 13 (2009). 

[17] A. W. Harrow and R. A. Low, Comm. Math. Phys. 291, 
257 (2009). 

[18] W. G. Brown and L. Viola, arXiv:0910.0913 (2009). 
[19] J. Emerson, Y. S. Weinstein, M. Saraceno, S. Lloyd, and 

D. G. Cory, Science 302, 2098 (2003). 
[20] A. J. Scott and C. M. Caves, J. Phys. A: Math. Gen. 36, 

9553 (2003). 



[21] L. Vandenberghe and S. Boyd, SIAM Review 38, 49 
(1996). 

[22] We used the the Matlab programs YALMIP and SeDuMi 

for solving the semidefinite programs. 
[23] I. Bengtsson and K. Zyczkowski, Geometry of Quantum 

State (Cambridge University Press, 2006). 
[24] S. T. Flammia, A. Silberfarb, and C. M. Caves, Found. 

Phys. 35, 1985 (2005). 
[25] J. Finkelstein, Phys. Rev. A 70, 052107 (2004). 
[26] S. L. Braunstein and C. M. Caves, Phys. Rev. Lett. 72, 

3439 (1994). 

[27] V. A. Osipov, H.-J. Sommers, and K. Zyczkowski, e-print 

arXiv:0909.5094 (2009). 
[28] S. Chaudhury, A. Smith, B. E. Anderson, S. Ghose, and 

P. S. Jessen, Nature 461, 768 (2009). 
[29] D. Gross, arXiv:0910.1879 (2009). 



